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1 Introduction 



Correlation functions are the main quantities characterizing quantum theories. From an experi- 
mental point of view, measurable physical quantities are directly related to dynamical correlation 
functions (or equivalently to their Fourier transform, the structure functions), for example in 
neutron scattering experiments on ferromagnetic crystals [H [2]. In the context of two dimen- 
sional integrable models various approaches to the computation of correlation functions and 
form factors have been developed over the years. It started first with models equivalent to 
free Fermions for which considerable works have been necessary to obtain full answer (see e.g. 
[31 m m m 13 El ) • Going beyond the free Fermion case has been a major challenge for the last 
twenty years. 

For integrable massive l-f-l-dimensional quantum field theories, form factors are accessible 
(using some hypothesis) from the bootstrap form factor program pioneered in the late 70's 
\10 \ [T2]: the analytical summation of their series corresponding to the correlation functions 
of local operators remains however an open question, although accurate numerical results exists 
for many theories and correlation functions. 

For integrable quantum spin chains [1'6 \ 1141 [TE\ and lattice models [16 1, the first attempts to 
go beyond free Fermion models relied on the Bethe ansatz techniques [17, and was undertaken 
by A. G. Izergin and V. E. Korepin (see e.g. [13] and references therein). Their approach yields 
formulae for the correlation functions [HI [El [20l [21] written as vacuum expectation values 
of some determinants depending on so-called "dual fields" which were introduced to overcome 
the huge combinatorial sums arising in particular from the action of local operators on Bethe 
states. However these formulae are not completely explicit, since these "dual fields" cannot be 
eliminated from the final result. 

In the last fifteen years, two main approaches to a more explicit computation of form factors 
and correlation functions have been developed, mainly for lattice models. 

One of these approaches was initiated by M. Jimbo, T. Miwa and their collaborators [221 
[2^ and enables, using some hypothesis, to compute form factors and correlation functions 
of quantum spin chains of infinite length (and in their massive regime) by expressing them in 
terms of traces of g-deformed vertex operators over an irreducible highest weight representation 
of the corresponding quantum affine algebra. These traces turn out to satisfy an axiomatic 
system of equations called q-deformed Knizhnik-Zamolodchikov (q-KZ) equations, the solutions 
of which can be expressed in terms of multiple integral formulae. Using these equations similar 
formulae can be obtained in the massless regime. Recently, a more algebraic representation for 
the solution of these q-deformed Knizhnik-Zamolodchikov equations have been obtained for the 
XXX and XXZ (and conjectured for the XYZ) spin 1/2 chains; in these representations, all 
elementary blocks of the correlation functions can be expressed in terms of some transcendental 
functions [25l [261 ET]. A detailed review of the approach can be found in [14J. 

The second approach has been developed by N. Kitanine, J. M. Maillet and V. Terras 
[28l [29l [30] . It combines the algebraic Bethe ansatz techniques [El [H] with the solution of 
the so-called quantum inverse scattering problem [2S\ I29j . It leads in particular to explicit 
determinant formulae for form factors of the finite Heisenberg spin 1/2 XXX and XXZ chains 
and to their correlation functions as explicit multiple sums. The solution of the inverse scattering 
problem means in practice finding an explicit realization of the local operators of a large variety 
of quantum spin chains in terms of the quantum monodromy matrix entries appearing in the 
algebraic Bethe ansatz framework and containing in particular the creation operators for Bethe 
eigenstates of the chain. Hence the computation of elementary blocks of the correlation functions 
reduces to a soluble algebraic problem in the Yang-Baxter algebra generated by the monodromy 
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matrix entries. Elementary blocks of correlation functions of the infinite XXZ spin 1/2 chain 
both in the massive and massless regime as well as in the presence of a magnetic field [30] have 
been computed in terms of multiple integrals. At zero magnetic field it gives a complete proof of 
the multiple integral representations obtained in [22l[2lj both for massive and massless regimes. 
Hence, together with the works [221124]. it also gives a proof that correlation functions of the XXZ 
(inhomogeneous) chain indeed satisfy (reduced) q-deformed Knizhnik-Zamolodchikov equations. 
In addition this method has proven to be effective in dealing with spin-spin correlation functions 
in the presence of a magnetic field [311 [5^ , dynamical correlation functions and at non zero 
temperature [35j, cases which were out of reach of the vertex operator method. A recent review 
of this approach can be found in [15] . 

So far, we have only mentioned various works dealing with spin 1/2 chains. However, consid- 
erable effort has been made to extend the results mentioned above to the higher spin integrable 
chains [351 \36\ [371 138] . Integral formulae for the correlation functions of the spin 1 XXZ chain 
have been obtained in [391 I^Ol E] • The correlations of spin chains of arbitrary spin have also 
been studied: integral formulae for the form factors of local operators of the XXZ chain were 
obtained in [32] and for the correlation functions of the XYZ chain in [l3]. From the Bethe 
ansatz method, the solution of the quantum inverse scattering problem has been given also for 
higher spin cases in [29] leading in [33] to integral formulae for the correlation functions of the 
spin 1 XXX spin chain. 

The aim of the present article is to employ the Bethe ansatz approach and the solution of 
the quantum inverse scattering problem for higher spins to obtain determinant formulae for the 
form factors of spin operators of general integrable Heisenberg quantum spin chains in arbitrary 
(mixed) finite dimensional spin representations. More precisely, we will consider integrable 
Heisenberg spin chains with different finite dimensional spin representations at each site. With 
regard to their physical properties (in particular the characteristics of the zero temperature 
ground state), these "mixed" spin chains can be subdivided into two groups: those where spins 
are mixed in fixed proportions and those where a particular type of spin dominates while other 
spins can be regarded as impurities. Examples of the first type of chains are alternating spin 
chains. The second type are spin chains with magnetic impurities. The eigenstates, spectrum 
and thermodynamic properties of both classes of models have been studied within the Bethe 
ansatz framework. The alternating spin 1- spin 1/2 chain has been studied in [451 1461 137]. More 
general combinations of spins have been dealt with in [481 149] . The spin 1/2 chain with one 
spin s impurity has been studied in [50] and especially in the context of the Kondo problem 
[511 152] . The more general situation of a spin s XXX chain with one spin s' impurity was 
analyzed in [53] (the special case s = 1 has been studied in [531 [55]). The effect of impurities 
in Heisenberg quantum spin chains has also been the object of experimental investigation, as 
shown for example in [56]. 

In this article we will compute the form factors of local spin operators for the general inte- 
grable XXX (mixed) spin chains and express our formulae in terms of determinants of elementary 
functions, similar to those found in [28] for the spin 1/2 case. 

This paper is organized as follows: In section 2 we review the algebraic Bethe ansatz frame- 
work for the general XXX quantum spin chains. In section 3 we derive closed formulae for all 
form factors of spin operators of the XXX chain in arbitrary spin representations. We employ 
our formulae to compute the total magnetization of the chain. In section 4 we give our con- 
clusions and discuss some perspectives. Some lengthy computations are given in appendices: 
appendix A contains proofs of two identities involving the higher spin eigenvalues of the transfer 
matrix which we have employed in our form factor computations and appendix B presents an 
alternative formula for the solution of the quantum inverse scattering problem for local spin 
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operators useful in form factor computations. 



2 Algebraic Bethe ansatz for XXX quantum spin chains 

In this paper we will consider general integrable XXX quantum spin chains of length N and 
periodic boundary conditions 5'^_,_j = S"^ for a = z, it. From the algebraic Bethe ansatz view 
point, the fundamental object characterizing such chains is the quantum monodromy matrix, 
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given as a tensor product of i?-matrices solutions of the Yang-Baxter equations |57[ [5H] which 
express the quantum integrability of the system. 



Rif'^\x 



(2.2) 



where V^^^ are (2s+l)-dimensional vector spaces, usually chosen as C^*"*"^. The quantities {^} = 
• • • 1 ^n} arbitrary inhomogeneity parameters attached to the sites of the chain. The spins 
si, . . . , Sat at each site of the chain are belonging to arbitrary finite dimensional representations in 
()2.ip which can be a priori different for different sites of the chain. The auxiliary space, indicated 
by the index, has been chosen to be 2-dimensional, which allows us to rewrite r(A) in the 
standard matrix form with operator entries {A(A), i?(A), C(A), I5(A)fl As a consequence of 
the Yang-Baxter equations for the i?-matrices, these operators satisfy a Yang-Baxter quadratic 
algebra and the traces of the above constructed monodromy matrices, namely the transfer 
matrices, 

t(i/2)(^) = (^ + ^)(^), (2.3) 



commute among themselves for arbitrary spectral parameters /U, A, namely. 



(2.4) 



leading to a large commutative subalgebra. Thus, the above constructed monodromy matrix 
describes a general XXX integrable quantum chain having different (finite dimensional) spin 
representations at each site of chain. 

A particular example of i?-matrices is the one corresponding to the spin 1/2 representation, 
the well-known 4x4 matrix 



where the functions 6(A) and c(A) are 



/ 1 \ 

6(A) c(A) 

c(A) 6(A) 

V 1 / 



(2.5) 



6(A) 



X + r] 



and c(A) 



X + r] 



(2.6) 



*We will usually not make explicit the dependence of the operators A, B, C and D and the monodromy matrix 
(|2.ip on the inhomogeneity parameters {^}. 
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and r] is an arbitrary parameter. The matrix (j2.5p is normalized in such a way that the following 
unitarity condition is automatically satisfied 



R 



12 l^J-n'21 



(-A) = I 



12, 



(2.7) 



where / is the identity matrix. i?-matrices and monodromy matrices of higher spin represen- 
tations can be constructed from (j2.6p by means of the fusion procedure developed in [35]. It 
gives, 

rjS+ 



R^^''\X) 



X + T]U-S 



(2. 



where S^,S^ are generators of the su{2) algebra in the spin s representation with standard 
commutation relations. 



2.1 The Bethe ansatz equations 

Having introduced R- and T-matrices we can now proceed to characterize the energy eigenstates 
of the chain. The key observation there is that the local Hamiltonian operator can be identified as 
a simple function of the transfer matrix. Hence finding its spectrum reduces to constructing the 
common spectrum of transfer matrices i^^/^) (/i) for arbitrary values of the spectral parameter n 
(since they all commute). The algebraic Bethe ansatz method starts by constructing a particular 
eigenstate of the transfer matrix, called the reference state |0), which is annihilated by the 
operator C(A). In order to identify such a state we notice that for every vector space Vj one can 
find a vector \0)j of dimension 2s j + 1 which fulfils the following equality 



1 



(2.9) 



Here, and in the following, we use the notation = u ±rj/2 for any complex number u. The 
reference state |0) is nothing but the tensor product 



N 

[0) = (g) \0)j with |0), 



/1\ 





(2.10) 



that is a completely ferromagnetic state with all spins up. The monodromy matrix (|2.1|) acts 
on the reference state as 

,(1/2) .,v,„, / 1 * 



namely 
with 



nr.(m = ^ ■ j 10). 

A{X)\0) = |0), C(A)|0) = and D{X)\0) = d{X)\0), 



j=i - s + j=i 



n biX - q -{k- Sj)r^) 



k=l 



(2.11) 
(2.12) 

(2.13) 



4 



It follows from the definition above that d{X) has zeros at A = '^^jV > ^i^d poles at A = — Sji] 
for all J = 1, . . . , A^. Let us now proceed to the description of the spectrum. Eigenvectors of the 
transfer matrix other than |0) will be constructed within the Bethe ansatz framework by acting 
successively on |0) with operators -B(Aj). In general such eigenvectors will be of the form 

\^{{X})) = B{Xi)...B{Xe)\0). (2.14) 

We want to select out those states (I2.14p which are eigenstates of the transfer matrix and hence 
of the Hamiltonian. The condition that (j2.14p is a common eigenstate of (|2.3p for all values of 
is a set of coupled algebraic equations for the spectral parameters Aj describing the state (I2.14p 
which are known as Bethe ansatz equations 

The eigenvalues of t^^^'^\p) corresponding to the eigenvectors ()2.14p fulfilling ()2.15p are given 
by 

t(V2)(^, {A}) = J] fe-i(A,. - ;.) + d(^) n - A,). (2.16) 

2.2 Fusion 

The mechanism of fusion for quantum spin chains was first developed in |35j for the XXX spin 
chain (and later on in [59] for the XXZ spin chain). It provides a procedure for constructing 
higher spin objects (i?-, monodromy and transfer matrices) starting with spin 1/2 objects. The 
fusion identities for the quantum monodromy matrix can be written as follows 

\^ * + (s - 1)??)7^(12) - r/) / 

where for the sake of simplicity we dropped the quantum indices 1, . . . , of the monodromy 
matrices ()2.ip and the monodromy matrix associated to the spin s in the auxiliary (first) space 
can be computed as a product over matrices as, 



7^o;liv(A;{e}) = 4lv'"^(A-e^)---<''^^(A-e,)---4r^(A-ei) 

with, 



(2.18) 



- io) ^ End(Fj^^ ® Vf'^). (2.19) 
The matrix 

is a direct sum of projectors onto the vector spaces 
resulting from the tensor decomposition 

^^^^^#-^(S©^(V' (2.21) 
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with Vr, ^ 



and V, 







Further 



X{u) = A{u+)Diu-) - B{u+)C{u-), 



(2.22) 



is the quantum determinant of the transfer matrix (|2.ip . which commutes with all operators 
A, B,C, D. A completely analog expression holds when replacing T-matrices by i?-matrices 
which permits the construction of (|2.8p starting with (j2.5p . A direct consequence of (j2.17p is 
the following recursive relation involving transfer matrices of various spins, 



and therefore, 



(2.23) 



(2.24) 



for the corresponding eigenvalues on a Bethe state |^({A})) and where for simplicity we have 
dropped the explicit dependency of the eigenvalues r on the spectral parameters A of the Bethe 
state. In [59] it was also realized that, as a consequence of ()2.23p . a generating functional for 
monodromy matrices of higher spin could be constructed as 



F{z) 



1 



oo 

^ zH^''/^\u- +kri/2). 



(2.25) 



Here we take t^^^x) = 1 for any values of x, and define z as a shift operator which acts on an 
operator t{x) as 

t{x)z = zt{x + 7]). (2.26) 

There are some differences between (I2.25P and the formula given in [59], which are due to the 
different normalization chosen for the i?-matrices in this paper. With the help of (|2.23p we can 
easily see that indeed ()2.25p generates any higher spin monodromy matrices. One must only 
match the terms on the l.h.s. and r.h.s. which contain the same powers of z, paying special 
attention to the fact that z is a shift operator. For example, expanding the l.h.s. of (I2.25P we 
find the following first terms 



F{z) = l + zt^'^/^'>{u) + 



zt 



(1/2)/ 



{u)zt^^l'^\u)-z^x{u^) 



+ ... 



(2.27) 



The 0(1) and 0{z) terms match trivially those on the r.h.s. of (|2.25p whereas the 0{z'^) term 



IS 



^i(l/2)(^)l _ ^2^(^+) ^ ^2 r(l/2)(^ ^ ^)t(l/2)(n) _ ;^(^,+ )l = zh^^\u+), 



(2.28) 



which exactly agrees with (j2.23p when taking s = 1. Similarly, one can establish the agreement 
with (I2.23P at every order in z. 

The existence of the generating function F{z) translates into the existence of a generating 
function for the eigenvalues of the operators t^^\u). They appear as coefficients of the corre- 
sponding Laurent series in z for the generating function. This allowed the authors of [59| to find 
explicit formulae for the eigenvalues which, translated into our present normalization, take the 
form 



2s 



r(^)(^,{A}) = j;c(^)(n)n 



{u^ - \p + sri){u ~Xp-sri) 



a=0 



p=i ^'^'^ - Ap + (a - s)rj){u - Ap + (a - s)ry) ' 



(2.29) 
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with 



2s-l 



Yl d{u+ + {k - s)r]) and C^J(u) = 1. 



(2.30) 



k=a 



It is now easy to check that (|2.29p indeed solves (|2.24|) . Similar relations hold for the XXZ case 
as well. 



2.3 The inverse scattering problem 



In the context of the algebraic Bethe ansatz approach, the inverse scattering problem is un- 
derstood as the problem of expressing quantum local operators of the chain in terms of the 
operators A,B,C,D with generate the Bethe states. The solution to this problem for a large 
variety of spin chains was provided in [28^1291. In particular, the reconstruction formulae for the 
spin operators of the XXX chain in an arbitrary spin Sj representation were found to be 



with 



.k=l 



2s 



tlt^'^H^k)-' 



.k=l 



a = ±, z 



(2.31) 



1) f n + ^) a3)(.- + ik - .),)t(^) (u- + ^-^^ 



k=i 

g t(^) (n- + i^i^) A^^in- + {k - .),)t(-l) + ^ 



(2.32) 
(2.33) 



The expression (j2.32p was derived in [29j using (j2.17p . whereas the non-trivial equivalence be- 
tween the expressions (I2.32p and (I2.33P is established in appendix B of the present manuscript. 
The existence of these two expressions will be very useful for later form factor computations. 



3 Form factors 

Let us start by introducing several formulae which we will use in the course of our computations. 
First of all, we will need the action of the operators A and D on a generic quantum state 
l^({A})) ■.= Ui=iB{Xk)\0), 



Aix)\^{{X})) = llb-Hx,-x)\^i{X})) 



k=l 



Yl c(Ap - X - r/) J] b-\Xk - \p)B{x) II BiXk) |0) , 

p=l kj^p k^p 



(3.1) 



Dix)\^{{X})) = dix)llb-\x-X,)\^{{X})) 



k=l 



+ Yd{Xp)c{Xp-x-ri)l[b-\Xp-Xk)B{x)llB{Xk)\0). (3.2) 



p=i 



k^p 



k^p 
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We see that this action is divided into two kinds of terms, which are traditionally refer to as 
"direct" and "indirect" terms. Direct terms are those which leave the original state unchanged 
up to a scalar factor while in indirect terms one of the parameters Ap has been replaced by the 
parameter x. 

Secondly, we need to compute scalar products of quantum states, either of two Bethe states 
or of one Bethe and one generic state. The formula for the scalar product of two Bethe states was 
originally obtained in [60ul9j . Later on it was proven [61j that a completely analog formula also 
holds for the scalar product of a Bethe state and an arbitrary state. Finally, the same formula 
was re-derived in [:28] with the help of the -F-basis introduced in [62]. It takes the following form, 



WM)I^({A})) :=5,(W,{A}) 



detg({/i},{A}) 

n(Ai - Aj)(/ij- -/ij 

i<j 



(3.3) 



where H({fj,},{X}) is a i x £ matrix of components 



H, 



ab 



V 



(3.4) 



and = {/ix, . . . , /U^} and {A} = {Ai,...,A^} are a Bethe state and an arbitrary state, 
respectively. 

Finally, we would like to recall a property of determinants which we will employ in the 
computation of the magnetization below: given two ^ x £ matrices X and y such that all rows 
of y are identical (rank 1 matrix) the following equality holds. 



+ y) = deiX + det X^^'^ , 



where 



"^ab 



Xab for bi^p, 

yap. 



(3.5) 



(3.6) 
(3.7) 



3.1 The form factors of S"! 

We define the non-vanishing form factors of the local operator in the spin Sj representa- 
tion as 



i^/(j,M,{A}) = (V(M)|5n^({A})) 



(3.8) 



with {/i} and {A} being two sets of d. Bethe roots, therefore characterizing two Bethe states. 

Inserting (I2.3ip with (I2.32p and a = z into (13. 8p we obtain the following sum of matrix 
elements 



^.0,M,{A})^^^E 

x(^(M)|(A-Z))(^,(A:))t(^) 



^7 + 



■ 2 
{k — 2sj)r] 



(3.9) 
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where the functions 



^■(M) = flr(^'=^(e.,M), 



(3.10) 



k=l 



originate from the actions 



k=l 



'A,-i(M)(^(M)l, 



(3.11) 
(3.12) 



k=l 



and we introduced the variable 

u,{k) = ^- + {k - sj)v. (3.13) 

In ()3.9p we have not yet replaced the operators by their eigenvalues on the Bethe state 

l^'jA}) for reasons which will become apparent below. Rewriting now A — D = 2 A — (A + D) 
and recalling the general action of the operator ^ on a Bethe state given in (|3.ip we find that 
()3.9p is equivalent to 



F/ j,{/x},{A} = ' 

24>A{X}) 



[5(j,{M})5KM,{A}) 



(3.14) 



2si e 



k=ip=i 



{7 + ^^^^,(A} 



V 



/ip - i^j(A:) 



k^p 



Here g{j, {//}) is the function defined in equation (jA.ip of appendix A. The term proportional to 
g{j, {^}) collects all contributions which are "direct" in the sense indicated after (I3.ip - ()3.2p . The 
form of this term results from acting with both operators t^^j "'^/^) and tC^"^)/^ on the Bethe state 
l^({^}))- By carrying out the computation in this way we achieve that the direct contribution 
to the form factor is proportional to the function g{j, {a*}) and we can take advantage of the 
identity (jA.lip proven in appendix A. 

The matrix element S'({A}, {^u, /ip — > Vj{k)}) in ()3.14p is the scalar product (13. 3p with the 
parameter replaced by Vj{k). It can be written as 



s{{ii,^i^ x},{\}) = II- \Trr\ — VY( 



(3.15) 



i<j 



where H^'p\x) is a matrix such that H^^^ = Hab{{X}, {/i}) ioi h ^ p and 



V 



Xa-X 

t{Xa,x) 
b{Xa - X) fj^ 



Y[iXi -x + r]) - d{x) W{Xi - x-rj) 



(3.16) 



n(A. 



U^-Hx. 



d{x)Y[b~^{x - Xi 
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with 



t{X,x) 



7] 



(3.17) 



{X-x){X-x + r]) 

and p(^/^)(A) being the momentum defined in (IA.6p . Inserting (j3.15p into (I3.14p with x = Vjik) 
we obtain the following expression for the form factors 



i^/(j,W,{A}) 



g{j, {/.}) det if - 2 E n (^fc - /^p + V) det HiP) (C,) 



p=l k=l 



20,({A}) 



■, (3.18) 



where 



H 



(p) 



Hab{{X},{l^}) for 
/(P)(j,{^},{A}) 



(3.19) 
(3.20) 



and f^^\j, {/i}, {A}) and g{j, {/u}) are the functions introduced at the beginning of the appendix 
A. Thanks to the identities (|A.3P and ()A.4p this expression can be simplified to 



i I 



F/(j,M,{A}) 



det ii - E n (/"fc - /^p + r/) det Z(f ) (C^.) 
'A,({A}) n(A.-A,)(^,.-^,) 

Kj 



(3.21) 



with 



ab 



^(P) 



i^a6({A},M) 
-idaP^^^\Xa-q) n 



for b ^ p, 

Afc - ^7 + Sj?? 



(3.22) 
(3.23) 



Notice that, since the entries zj^ depend on both indices a,p we can not write (|3.2ip in terms 
of a single determinant by employing ()3.5p . However, the factor 111=1 (/^fc ~ ^p + ^) '^i^ be 
canceled when normalizing by the norm of the Bethe states, which allows a particularly simple 
form for the magnetization. We compute this quantity in the next section. Let us also remark 
that the computation of (|3.8p can be done in two different ways, which of course must lead to 
the same result: in ()3.14p we decided to make the operators {A — D){uj{k)) act on the Bethe 
state however we could as well have chosen to act on the Bethe state {A}. If we do that 
and employ the same reconstruction formula ()2.32p as before, we will obtain an expression which 
looks rather different from (j3.2ip and cannot be simplified in any obvious way. However, we 
can alternatively employ the reconstruction formula p.33p to do the same computation. In that 
case we will obtain a formula completely analogous to p.2ip where the roles of {A} and {/x} are 
just exchanged. 



3.2 Magnetization 

One interesting physical quantity we can now compute is the magnetization at site j which is 
defined as 
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The norm of the Bethe state is given in terms of the Gaudin matrix $'({A}) as fohows |60j 

(3.25) 



mWmW)) = 7?^n^"HAa-A,)detcI>'({A}), 



with 



( 



1 ^ 



b{\a - Afc) 
d{\a) J-^A h{\k - \a) 



\ 



(3.26) 



and can be obtained from the scalar product formula (|3.3p in the limit \a — > /i^ for all a = 1 . . . £. 
In particular, it is easy to prove that for two identical Bethe states, the entries of the matrix 
H{{X}) defined in (13. 4p become proportional to those of the Gaudin matrix as 



Hab = -\b + r]Wab- 



and therefore 

det H{{\]) = \[{\- Xj + ri) det $'({A}). 
Using these identities and (|3.18|) we obtain 



(3.27) 



(3.28) 



(5| 



where the matrix A^(") is such that 



Sjdet $'-Ea=i det 



det 



M 
M 



(a) 
(a) 



^ij for i / a, 



(3.29) 

(3.30) 
(3.31) 



where p^^^\X) is the momentum (jA.6p defined in appendix A. Since A^j" is independent of the 
value of the index i, the identity (13. 5p can be used to bring (I3.29P into the form 

det{^' -M/sj) 



det 



(3.32) 



where M is a rank 1 matrix whose entries are Mab = • The total magnetization of the 
chain can be computed as 



N 



where Al+ot is a rank 1 matrix with entries 



det{<^' - Mtot/sp) 
det 



(3.33) 



TV 



(Mot)ia = 5]-^9A,p(^^)(A,-e-), 



(3.34) 
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and So is the total spin of the completely ferromagnetic state |0) 



N 

so = ^Sj. (3.35) 

i=i 

The total magnetization can be computed by noticing the following property of the Gaudin 
matrix (fXM]) 

e 

^<, = (Mot)fca. (3.36) 
6=1 

Therefore 

Pttot = ■So det (^I = So- I, (3.37) 

where / denotes the identity matrix and U is the matrix, 

Uab = l V a, b. (3.38) 

This computation provides a consistency check of the general formula ()3.18p as it reproduces 
the values of the magnetization known for particular cases. For example, for the pure spin 1/2 
XXX chain in the ground state we have £ = N/2 and sq = N/2 so that the total magnetization 
vanishes (see e.g. )28j). Another well known example is the spin 1/2 XXX chain with one spin 
s impurity. In this case the ground state is also characterized by ^ = N/2 roots, but the total 
spin of the chain is sq = (-^ ~ l)/2 + s. Therefore the total magnetization becomes s — 1/2, in 
agreement with the value computed in |52j . 

3.3 The form factors of 5+ 

In this section we compute the form factors of the operator Sj' of the quantum XXX spin chain. 
That is a spin generator sitting at site j of the chain and living in the spin sj representation. 
The only non- vanishing form factors are 

where |^'({;u})) and |^'({A})) are two Bethe states with i + 1 and i Bethe roots, respectively. 
Prom (f23T]l and (|232]) we obtain 

F^U, {A}, {,}) = ^f0^ E b'''-^' fe + ^W2, {A}) 

xt(^) (ej + (A; - 2sj)r,/2, {^}) 5,+i({/i}, {A, Uj{k)})\ , (3.40) 

where I'jik) is again the variable introduced in (|3.13|) and S'^+id/x}, {A, z^j(/c)}) is the scalar 
product introduced in ()3.3p with {A, i^j{k)} the set of ^ + 1 variables {Ai, . . . , A^, i^j{k)} in the 
given order, that is I'jik) is the variable number i + 1 (notice that the ordering of variables is 
important in ()3.3p ). Employing the definition ()3.3p we have 

S,,M^}.{X.^m) - n(A. - (3,41) 

n (Mj - M n (Ai - Aj) fj- 

l<i<j<i+l l<i<j<^ 
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where H is a matrix with components Hab given by (j3.4p . except for the column £ + 1 whose 
entries are given by (13.160 with x = vjlk). This means that ()3.40p can be rewritten as 

(hi i(|A|) deti?+(£,) 
<t>j{W) 11 11 (Ai-Aj) 



where 



J/+ = {A}) for + (3.43) 

^i+i = /^'^^(j,{A},{m}), (3.44) 

with f^"'\j, {A}, {/u}) being the function (|A.2p defined in the appendix. Exploiting the identities 
proven in appendix A we can bring the form factors into the simple form 

£+1 

r,+rs^^s^^ '^.-i(W) "i^^' " " '^'"^^ detC(g,) 

FiO^ {A}, M = / .. —f ff 7 — 7T r^, 3.45 

0j_i({m}) /f ^_ . n n (A»-Aj) 

11 K'^k — — Sjl]) l<i<j<i+l i<«<i<^ 

fc=l 



with 



Cab = i/afe(W,{A}) for b^£ + l, (3.46) 
Ca^+i = -i5^y^^)(^,-e7). (3.47) 

3.4 The form factors of Sj' 

In this section we compute the form factors of the operator of the XXX spin chain in the 
spin Sj representation 

F,-{j,{f^},W) = m{f^})\ Sr \^ax})) . (3.48) 

We can compute these form factors along the same lines of the previous section. However if we 
take (j2.3ip with (j2.32p as starting point for our computations we will obtain closed formulae for 
(I3.48P which lack the simplicity of (I3.45p . It is in fact more convenient to use (I2.33P instead of 
()2.32p . Using formula (j2.33p and proceeding as in the previous section it is not difficult to prove 
that the form factors (|3.48p are related to (|3.45p as 

4 Conclusions and outlook 

In this paper we have obtained general expressions for the form factors of all spin operators 
{S'^,S'^} of the integrable XXX quantum (higher) spin chain. Our formulae hold for any spin 
representation of the operators as well as for any spin configuration at the remaining sites of the 
chain. These results can be extended to the anisotropic XXZ spin chain as well. 

In view of the very general nature of our formulae we expect they will be useful for the study 
of a number of interesting physical systems, whose thermodynamic properties have been already 
extensively studied in the literature but whose correlation functions remain unknown. Amongst 
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these systems alternating spin chains, first studied in [ISldHlllH] and more recently in [15], and 
impurity systems, such as the Kondo model [5ll [52] and the systems considered in [50 [ [53l [63] . 
seem specially interesting examples. 

Recent results for the XXX and XXZ spin 1/2 chain show that determinant formulae for 
the form factors, similar to those obtained here, can be successfully employed for numerical 
computations of spin-spin dynamical correlation functions of finite chains, which can match 
very precisely experimental data [64> .65j . The form factors computed here might serve for the 
numerical study of the correlation functions of impurity systems. 

From the analytical point of view, the most natural continuation of this work is the com- 
putation of correlation functions for mixed spin chains along the lines of [3^ [HTl [32] • Further 
interesting generalizations of this would be the study of dynamical correlation functions in the 
spirit of [33] and the computation of correlation functions at finite temperature, generalizing the 
program initiated in [34j . 
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A Proof of two identities for the eigenvalues of higher spin 
transfer matrices 

In this appendix we prove two identities allowing for a compact formula for form factors. Let us 
consider two Bethe states {A} = {Ai, . . . , A^} and {/i} = {/x^, . . . , /x^} with t and i Bethe roots 
respectively and define the functions 



2s, 



5(i,{A}) = ^ 



fc=l 



p=l p=l 



(Al) 



and 



2s,- 



k=l 



/('^)(i,{A},M) = ^ 



{k — 2sj)r] 



{/.} (A.2) 



n b-\^,^ - uj{k)) - d{u,{k)) n b-\u,{k) - /.p) 

p^a p^a 



n {^^p - ^ji^)) 



p=l 



\[{\p-Vj{k)) 
p=l 
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with I'j^k) given in (j3.13p . The aim of this appendix is to prove that 



with, 



and 



;(i,{A}) = 2.,r(^^)(6,{A}), 



/('^)(j,{A},M) 



n (Pp - ?j + sjv) 2s, 

n(Ap-^7-s,r?)^'=i 
P=i 



2s, 



k=l 



A — SjT] 



p(«^)(A) =ilog 



A + Sji] 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



is the momentum of a spin sj pseudo-particle of "rapidity" A. In particular, the identities above 
imply 



/W(j,{A},{A}) =r(^^)(C,,{A}) \-id,J-^^\X, - C) 



(A.7) 



Let us commence our proof by writing down the explicit expression of the eigenvalues following 
from eq.^M) 



2sj k 



n 



^ a=0 
jS.^ + SjT] - Xp){iyj{k) - Xp) 

L {uj{k) + {a + l)r] - Xp){uj{k) + arj - Xp) 



{k — 2sj)r] 



k-l 



(3=0 



(i^) / {k-2s,)v 



n 



{Uj{k) - llp){ij - SjT]- lip) 



(A. 



(A.9) 



whose product enters the two formulae we want to prove. Before we insert these formulae in 
(|A.ip - (|A.2p it is useful to note that in (|A.8p the only non vanishing term is the one corresponding 
to a = 2s j — k. It is easy to see that any other terms will contain a factor d{S,J + Sjij) = and 
never contain the singular factor ~ ^jV)- Likewise it is straightforward to argue that such 
factors will never appear in cj^^^~^^^^\Cj + ik-2sj)v/2) for any of the allowed values of /3 and 
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k, and therefore each term in the /3-sum is non- vanishing and non-singular. Therefore 



^(^.-1) . + {A}) r(^) Uj + {k - 25,)r//2, {^} 



p=i 



k-l 



Y^c'^p"' \q + {k-2s,)ii/2) 



n (/^p - T^j{k)) 
p=i 



p=i 



(A.IO) 



n if^p - e+ + (sj - p)v){f^p - (J + (sj - m 
p=i 



where we used the second property stated in (j2.3Up . i.e. C2g^ (^^ + k'q/2) = 1. Using these 
formulae, the equalities (|A.ip and (|A.2p are equivalent to proving : 



EE 

k=l f3=0 



(,Xp-ujik)y^ 



2 J ^J,H\p-^^ + {sj-P)r,) 



n ^"'(Ap - - d{v^{k)) n - \p) 

p=i p=i 



2s,, 



(A.ll) 



and 



2sj fe-1 

EE 

fc=l/3=0 



(fc — 2sj)r] 



{^l„-Vj{k)f 



^,'h{^^p-^-t + {s,-m 



n - v,{k)) - d{v,{k)) W b-\uj{k) - f,^) 

_p^a p+a 



(A.12) 



respectively, with 



h{\):={\){\ + r]). (A.13) 

We commence by establishing (jA.lip . Let us consider those terms in (jA.lip which do not contain 
any factors proportional to the function d{x). The only terms of this kind correspond to taking 
(3 = k — 1 and only the first contribution on the second line of (|A.lip . that is 



2s, I 

En 

k=l p= 



(Ap-e, +{s,-k)r]){\p-r 



Sj - l)r]) 



2s 



J' 



(A.14) 



16 



therefore it remains to prove that ah remaining terms cancel each other. Those remaining terms 
are 



where we have used the identity 

Replacing the sum over A: by a sum over A;' = — 1 in the first line of (1A.15|1 we see that both 
terms are essentially identical but for the extra term k = 2s j contributing in the second sum. 

(s •) _ 

However such term is vanishing since Cy'{i~ + {k- 2sj)r]/2) = and therefore we have proven 
that indeed all terms containing the function d{x) in (jA.lip cancel each other. With this we 
have established (lA.lip and therefore (lA.ip . The proof of the identity (|A.12p can be carried out 
by following exactly the same steps, namely considering separately those terms which contain 
(i- functions and those which do not. 



B An alternative formula for the solution of the inverse scat- 
tering problem for XXX spin chains 

The solution of the inverse scattering problem for XXX spin chains and arbitrary spin represen- 
tations was given in [29]. This solution was recalled in formulae (j2.3ip and (|2.32p and employed 
in order to obtain closed expressions for higher spin form factors of the operators and S~^. 
However we have observed that in some cases the form factor formulae obtained from (12.320 
simplify considerably (thanks to the identities proven in appendix A), whereas in other cases 
such simplifications do not occur. For example, a simple formula for all form factors of the 
operator can be obtained form (|2.3ip - (j2.32p . but the formula we would obtain for the form 
factors of S~ from ()2.32p is much more complicated. This asymmetry seems rather unnatural 
and one may suppose that the formulae for the S~ form factors can be further simplified and 
recasted into analogous formulae as those for the form factors. This should be doable by 
using properties of determinants as well as the Bethe ansatz equations; however it turns out to 
be a rather non-trivial proof. An alternative way of finding simpler formulae for the form factors 
of S~ is to start with a different (but equivalent) reconstruction formula, that is (j2.33p . In this 
appendix we will show that (|2.32p and (|2.33p are indeed equivalent. The key tool are once again 
the fusion relations ()2.17p and one particular property of these noticed in ^35j . 



B.l An alternative version of the fusion relations 



Let us start by decomposing ()2.17p into the following two equations 

K.iayTl^h-- + sv)Tl:y'^\x-)P:^^^^ = t£\„^^(x), (B.l) 
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With respect to (j2.17p we have sUghtly changed our notation and introduced the index {a} to 
indicate that the space vj^^^ ^^'^^ ~ C^'^""'^ is isomorphic to the space V^^'*^ ^2^^) resulting from the 
fusion of 2s — 1 spin 1/2 quantum spaces 

vlpr^C^ j = 2,...,2s. (B.3) 



We can now employ fusion successively in order to express (|B.ip solely in terms of spin 1/2 
quantum monodromy matrices. By doing so, we obtain the following expression 



2s 



K.-.a^.U^ah^-' + (S - j)v)PL.a,. = ^^^..^^^ (x) . (B.4) 



It was proven in [35] that, for XXX spin chains, the operator Paj^^^^a2s ^ complete symmetrizer 
on all indices ai,...,02s. Let us now consider the formula (lB.2p . We can also express this 
relation entirely in terms of spin 1/2 monodromy matrices 

2s ^ 

PaAa}PL.a,. H + " J + l)^)^at ...a.. ^'.^ ^ = X{x + «r/)7;^\) (x) . (B.5) 

i=i 

As before the operator Pa2...a2s completely symmetric on the indices {02 • • • fl2s}- On the other 
hand the operator ^'"i^i is antisymmetric on the index ai (see |35J). The identities (|B.4p - ()B.5p 
can be further manipulated by employing the following property of the projectors which was 
stated in [35] 

This property allowed the authors of to prove that 

PL.a2.Il^af\x^ + {S - j)v)PL.a2s = P^.a^sU^^^ " " J»^a+ ...a.. , (B.7) 

and 

2s 

Pai{a}Pa2...a2s JI ^ij^^^ + ~ j)'n)Pa2...a2sPai{a} 
2s 

= Pa.ia}Pa2...a2.Il^^y'Hx- " " j)^)^'at ...a., ^''m • (B.8) 

The equalities (|B.7p - ()B.8p imply that the fusion relation (|2.17p still holds if in the arguments on 
the r.h.s. t] is replaced by —t], that is 

Pi2t[^\-^ - sv)Tt^'\x^)Pi2 = ( , ^ , ,,\^is-i). , ] ■ (B.9) 

y * X{x + {s - l)r))T^^2-^ ' {x - r]) J 

This equivalence was stated in equation (18) of [35] for i?-matrices, but it is easy to prove that 
it implies a similar property of the monodromy matrices. 
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The s = 1 case: As an example, let us prove (|B.7p for s = 1. In this case ()B.4p becomes 
simply 

P^,a,Tihx^)Tihx-)P^,a, = (B.IO) 

and the relation (|B.6p can be written as 

RiTal\±nv)=±P^,a,, (B.ll) 

for n = 1. The iJTT-relations 

Rifai\x)Ti}\x + /i)rij^(^) = ril V)rif^(A + ;u)i?iM^(A), (B.12) 
at A = 7/ and ji = x~ imply therefore 

Pla,Ti\x^)Tlk\x-) = ri|\x-)rif\x+)P+,,. (B.13) 
Multiplying (jB.lSp from the right and from the left by Paia2 obtain 

^'ia,rif^(x+)rij)(x-)p+,, = P+,,ri|^(x-)rif (b.u) 

which is precisely (IB.7P for s = 1. 

B.2 An alternative formula for the solution of the inverse problem 

In order to solve the inverse problem we need to find a systematic way of computing the traces 

Ai^^^ := Tro {s^T^il Miij)) , a = ±,z. (B.15) 

To do that, the key idea [29j is to use the fusion relations ()B.9p in order to obtain recursive 
formulae similar to ()2.23p for (IB.ISP . Recall that for the general XXX chain in finite dimensional 
representations, the co-product for local spin operator is trivial and therefore we can write the 
spin matrices in the two following ways 

= /i + S'^ ^ /2 ~ 5fi2) © 5fi2), a = ±,z. (B.16) 

Multiplying (|B.9p by S'fj and taking thereafter the trace over the spaces 1 and 2 on the l.h.s. 
and over the fused spaces (12) and (12) on the r.h.s. we obtain the following relations 

= A^^\x) + x{x + is- l)r?)A(f-i)(x - tj). (B.17) 

In addition, taking directly the trace on ()B.9P we obtain also new fusion relations for the traces 
of the transfer matrices 

t(m{x+ - si])t'^'-^/^\x+) = &\x) + x{x + {s- l)r])t^'-^\x - 7]). (B.18) 

The recursive relations ()B.17p are solved by 
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This can be proven by substituting Aq and in (|B.17p by the corresponding formulae 

from (1B.19|1 and using thereafter the fusion relations (jB.lSp . This is exactly the same solution 
(|2.32p up to the replacement rj — —rj in the arguments. By introducing a new index p = 2s — k+l 
we can rewrite (IB.19P as 

AL^)(n) = f; t(^) (u- + i^^) AS^(n- + (p - s)rj)&-'^^ (. + f ) , (B.20) 
p=i ^ ^ 

which is identical to p.32p up to the exchange of the transfer matrices. 
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